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I. INTRODUCTION 



Symmetry occupies a central and indispensable role in 
physics. Amongst others, it is the key concept behind the 
theory of phase transitions and critical phenomena. Most of 
the known phase transitions break one or more symmetries. 
When symmetry breaking occurs, the mean value of a quantity 
which is not invariant under the symmetry group of the theory 
becomes "spontaneously" non-zero in the broken symmetry 
phase. Such a quantity is called the order parameter. As well 
appreciated, the dimension D of the system and the symme- 
tries of the order parameter characterize different universality 
classes and play a pivotal role in determining whether a partic- 
ular phase transition may or may not exist. In particular, zero 
or one dimensional theories with short range interactions can- 
not exhibit a phase transition at any finite temperature. Addi- 
tionally, the Mermin- Wagner theorem [ 1 ] states that a contin- 
uous symmetry cannot be spontaneously broken at any finite 
temperature for two-dimensional theories with finite range in- 
teractions. On the other hand, Elitzur (2l demonstrated that 
a spontaneous breaking of a local symmetry is not possible. 
Below, we will show that Elitzur's theorem is a consequence 
of a reduction to zero of the effective dimension of the gauge 
invariant theory. Moreover, we will show that from the point 
of view of the non-invariant gauge fields, the presence of a "d- 
dimensional gauge symmetry" (see definition below) reduces 
the effective dimension of the theory from D to d. The dimen- 
sion d is intermediate between local gauge systems (d = 0) 
and global gauge symmetries (d = D). 

How is it possible that symmetry not only determines the 
nature of the transition but also its effective dimension? The 
answer to this question requires a definition of what we mean 
by local and "<i-dimensional gauge symmetries". A local 
gauge symmetry has its roots in a group of operations which 
leave the theory invariant and in which the minimum non- 
empty set of fields influenced by any of these operations oc- 
cupies a finite spatial region. A "d-dimensional gauge sym- 
metry" is any symmetry operation which leaves the Hamilto- 
nian (or action) invariant such that the minimum non-empty 
set of fields which are changed under the symmetry opera- 
tion occupies a d-dimensional region. In accord with this def- 
inition, a local gauge symmetry corresponds to the particular 
case d = 0. The above question is motivated by the recent 
interest in a considerable number of theories that are invari- 
ant under a d-dimensional gauge group with dimension d > 1 



S 11 II B 13 H 13 El As was already noticed SHU, 
some aspects of these theories indicate that the dimension is 
effectively reduced. In this paper we formalize this notion by 
proving a generalized Elitzur's theorem for any value of d. 

We will prove that the mean value of a local quantity which 
is non-invariant under a d-dimensional gauge group can be 
bounded by another mean value of the same quantity which 
is computed for an effective d-dimensional theory. The gauge 
transformations become a global symmetry of the effective 
theory and the range of the interactions in the original theory 
is preserved. This proof formalizes the idea of "dimensional 
reduction" and gives rise to four important corollaries. The 
first corollary is Elitzur's theorem. The second states that a 
d = 1 gauge symmetry cannot be broken at any finite temper- 
ature if the the theory only includes finite range interactions. 
The third corollary is a natural consequence of the Mermin- 
Wagner theorem: a continuous d = 2 gauge symmetry can- 
not be broken at any finite temperature if the interactions are 
of finite range. The fourth corollary extends Elitzur's theo- 
rem to certain systems in which at least an emergent low en- 
ergy gauge-like symmetry appears. By emergent symmetry 
we mean a group transformations which become a symme- 
of the Hamiltonian only within its lowest-energy subspace 
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Recently, theories containing d-dimensional gauge symme- 
tries with d > 1 became the focus of attention in different 
areas of physics. The Kugel-Khomskii model that describes 
the interplay between spin and orbital degrees of freedom in 
transition metals as well as orbital only variants embodying 
the Jahn-Teller interactions often display d = 1, 2 gauge-like 
symmetries |@, IE [9[ fToll . The well known planar orbital 
compass model laBGlJ] is invariant under a d = 1 Z2 -gauge 
group. The two-dimensional frustrated Ising magnet studied 
by Xu and Moore ^ provides another example of a d = 1 
Z2 -gauge invariant theory. There are also examples of frus- 
trated magnets fl2ll for which the d = 1 Z2 -gauge group is 
an emergent symmetry lITTIl . Here, the dimensional reduction 
may occur at low energies. Similar symmetries appear in ring 
exchan ge B ose metals fl3ll . Quantum Hall liquid crystalline 
phases iQ, fl5ll . lipid bilayers with intercalated DNA strands 
lll6ll . and "sliding Luttinger liquids" 0. 

The physical roots of our theorem and its corollaries are 
precisely the same fluctuations which generically inhibit order 
in low dimensions. Just as domain wall entropy eradicates any 
viable order in one dimensional systems having short range in- 
teractions, domain walls on different sliding phases of many 
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of the above examples inhibit discrete d=l gauge symmetry 
breaking at any finite temperature. Similar considerations ap- 
ply to continuous d= 1,2 dimensional gauge symmetry sys- 
tems whose order is destroyed by low energy excitations. 



II. A GENERALIZATION OF ELITZUR'S THEOREM 

We will start by considering the case of classical theories. 
As we will see later, the generalization to the quantum case is 
very simple. Before enunciating our theorem it is convenient 
to give a more precise formal definition of a d-dimensional 
gauge group. From now on we will consider theories which 
are defined on a lattice A. The extension to the continuum 
case is straightforward. A d-dimensional gauge symmetry of 
a theory characterized by a Hamiltonian H is a group of sym- 
metry transformations of H such that the minimal non-empty 
set of fields that is changed under the group operations oc- 
cupies a d-dimensional subset C C A. The index i denotes the 
sites of the lattice A. For instance, if a spin theory is invariant 
under flipping each individual spin the corresponding gauge 
symmetry is zero-dimensional or local. Of course flipping a 
chain of spins is also a symmetry, but the chain is not the min- 
imal non-trivial subset of spins that can be flipped. In general, 
these transformations can be expressed as: 



XJlk = Y[ gifej 



(1) 



where C\ denotes the subregion I, C\ C A, and A = (J z Ci- 

Theorem: The absolute mean value of any local quantity 
(i.e. involving only a finite number of fields) which is not in- 
variant under a d-dimensional gauge symmetry group G of 
the D-dimensional Hamiltonian H is equal or smaller than 
the absolute mean value of the same quantity computed for a 
d-dimensional Hamiltonian H which is globally invariant un- 
der the group G and preserves the range of the interactions in 
H. Non invariant here means that the quantity under consid- 
eration, /(0i) , has no invariant component: 



= o. 



(2) 



For a continuous group Eq.0 has to be replaced by 
J /[gi(0i)]^g = 0- To determine if there is a spontaneous 
symmetry breaking in the thermodynamic limit, the mean 
value of f(<fii) has to be computed in the following way: 



(ZOO) =lim h ^olim N ^ 00 (f(<l>i))h 1 N, 



(3) 



where (f(<f>i))h,N is the mean value of f{<j>i) computed on 
finite lattice of N sites and in the presence of a symmetry 
breaking field h. Since A = \J t Ci, the site i belongs at least 
to one set Cj . It is convenient to rename the fields in the fol- 
lowing way: <j>\ = ipi if i ^ Cj and <j>\ = rji if i G Cj.The mean 



value (f (fa)) h,v is given by: 



(f(<k))h,N 



E 



with, 



(4) 



(5) 



From Eq.®: 



\(f(fr))h,N\<\ 



I, (6) 



where {{[)} is the particular configuration of fields that max- 
imizes the expression between brackets in Eq.Q. H($,rj) is 
a d-dimensional Hamiltonian for the field variables 77 which is 
invariant under the global symmetry group Gj of transforma- 
tions XJjk over the field 77. We can define S(rj) = H(ip, 77). 
The range of the interactions between the ?7-fields in H(tj) is 
clearly the same as the range of the interactions between the 
^-fields in H{4>). This completes the demonstration of our 
theorem. Note that the "frozen" variables ipi act like external 
fields in B(rj) which do not break the global symmetry group 
of transformations XJjk- 

Corollary I: Elitzur's theorem. Any local quantity (i.e. in- 
volving only a finite number of fields) which is not gauge in- 
variant under a local or zero-dimensional gauge group has a 
vanishing mean value at any finite temperature. This is a di- 
rect consequence of Eq.Q and the fact that H (77) is a zero- 
dimensional Hamiltonian. 

Corollary II. Any local quantity which is not gauge in- 
variant under a one-dimensional gauge group has a vanishing 
mean value at any finite temperature for systems with finite 
range interactions. This is a consequence of Eq.© and the 
absence of phase transitions for one-dimensional Hamiltoni- 
ans with finite range interactions. Note that H (77) = iJ 77) 
is a one of these Hamiltonians and f(rji) is a non-invariant 
quan tity under the global symmetry group Gj [see Eq.©], 

Corollary III. Any local quantity which is not gauge invari- 
ant under a two-dimensional continuous gauge group has a 
vanishing mean value at any finite temperature for systems 
with finite range interactions. This results from the combina- 
tion of our theorem [Eq.©] with the Mermin- Wagner theo- 
rem: 



lim h ^lim N — v/lj " x ' 7 0. 

(7) 

Where we have used that Gj is a continuous symmetry group 
of H(rj) = Hty, 77), f(rj{) is a non-invariant quantity for Gj 
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[see Eq.©], and H(rj) is a two-dimensional Hamiltonian that 
only contains finite range interactions. fl9ll 

The generalization of our theorem to the quantum case is 
very straightforward if we choose a basis of eigenvectors of 
the local operators that are linearly coupled to the symmetry 
breaking field h. In this basis, the states can be written as a 
direct product \<p) = \i/j) \rj). Eq.© is re-obtained with the 
sums replaced by traces over the states \rj): 



IWi)>Mv|< 



Tr {rii}f( 7 H) e 



TT im} e 



-/3if(^,f?)-/3/iEi e c.*li 



(8) 



In this case, \tjj) corresponds to one particular state of the ba- 
sis \ip) that maximizes the right side of Eq.©. Generalizing 
standard proofs, e.g. I20I1 . we find in the quantum arena a zero 
temperature extension of Corollary III whenever there is a gap 
in the excitation spectrum: 

Corollary IV. If a gap exists in a system possessing a d < 2 
dimensional continuous symmetry (or, at least, an emergent 
continuous symmetry) in its low energy sector, then the ex- 
pectation value of any local quantity which is not invariant 
under this symmetry, strictly vanishes at zero temperature. 



III. PHYSICAL REALIZATIONS AND IMPLICATIONS 

Armed with the above results, we briefly review known 
consequences of Elitzur's theorem and then focus at length 
on several of the applications of our new results. 



A. Local gauge symmetries- Pure and Matter coupled gauge 
theories and Spin glass systems 

a) Gauge theories- The application of Elitzur's theorem in 
gauge theories is well known and reviewed in excellent works 
such as I2II1 . No work on Elitzur's theorem is complete with- 
out a mention of its most prominent application. In theories 
of matter at lattice sites (c^) coupled to gauge fields (Uij) re- 
siding on links, the action is a sum of a plaquette product of 
the gauge fields and (minimal) coupling between matter fields 
through the gauge 1231 . To illustrate, consider the action in a 
Z2 setting, 



S = —K^^UijUjkUkiUki - J y^^o-jUijcrj, 



(9) 



with Uij = ±1 and a, 



±1 elements of Z2. Many pioneer- 
ing results were found by l23ll . The action is invariant under 
the local gauge transformation G{ — > rji(Ji,Uij — > rjiUijTjj 
with rji = ±1 an element of Z2. Defining the gauge in- 



variant link variables Zj 



cfiUijCFj of which the action is 



a functional, we note that any correlator involving a prod- 
uct of any number of z's ("mesons") need not vanish. In 
the absence of coupling to matter (J = 0)- the "pure gauge 
theory"- only products of Uij along closed "Wilson" loops 
(W = (n*jec Uij)) are gauge invariant and may attain finite 
expectation values. 



b) Spin glasses- Perhaps one of the simplest realizations of 
spin-glasses is the Edwards-Anderson (EA) model f24ll with 



H 



^ ] J ij Si ■ Sj h • ^ ^ Si 



(10) 



Here, i and j are nearest neighbor sites of a regular lattice and 
the three component spins S = (S Xj S y ,S z ) are considered 
as classical vectors. In the EA model, the exchange constants 
Jij are independent random variables drawn from a (Gaus- 
sian) probability distribution P(J^ ). In spin glass systems, 
the experimentally measured expectation value of any observ- 
able corresponds to a quenched average over the distribution 



{ f {{S})) quench 



J dJP(J)(f({S}))j. (11) 



Here, (/({£})) j denotes the expectation value of f({S}) for 
a given quenched distribution of exchange constants {Jij}. 
The basic premise of spin glass systems is that the free en- 
ergies (and their derivatives) are to be averaged via the dis- 
tribution P(J) to obtain the corresponding quenched quanti- 
ties. The quenched expectation value of Eq.JTTl is, for any 
even distribution P(J), invariant under the following local Z2 
gauge transformation, 

Si ► TJiSi, Jij ► TjiJijTjj^ (12) 

with rji = ±1 (an element of Z2). 12^1 Performing the 
quenched disorder average of Eq.JTTl over J^ is a central task 
in spin glass problems and ingenious schemes have been de- 
vised. Parisi showed that the quenched average may be done 
by introducing n (with n — > 0) replicas of the spin field 
{Si}2=i at each site all coupling via the same J^ thus still 
satisfying the gauge transformation of Eq.(fT2b detailed above 
with a replica independent rji 1261 . Within the replica formu- 
lation, the average over J^ may be performed first leading to a 
functional of the replicated spin fields alone. To attain a finite 
expectation value, the quantity of interest must be gauge in- 
variant. The lowest order quantity invariant under local trans- 
formation is |3, 13 Qf = (Sf ■ S%) (as both spin field 
suffer a factor of r\i under the local gauge transformation, and 
rjf = 1). In this manner, we see how the usual spin-glass 
overlap order parameters (Q^ b and their sums) are dictated by 
Elitzur's theorem for the n-replica action, c) Many other ex- 
amples are found in numerous spin systems of current inter- 
est, For instance, Kitaev's 12911 exact solution of a special spin 
system on a honeycomb lattice owes its existence to gauge- 
like local symmetries implemented via spin products around 
individual hexagons. When these local symmetries are fused 
with Elitzur's theorem, we find that on-site magnetization is 
precluded. 



B. Exact Higher-Dimensional Gauge-like symmetries: Orbital 
systems 

a) Spins- In 3d orbital systems on cubic lattices, each tran- 
sition metal atom is surrounded by an octahedral cage of oxy- 
gens. Crystal field splitting lifts the degeneracy of the five 3d 
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orbitals of the transition metal to that of two higher energy 
e g levels (residing in the space spanned by the two orbital 
states \d 3z 2_ r 2) and \d x 2_ y 2) wherein the electronic wave- 
functions point towards the surrounding oxygens with an as- 
sociated high energy Coulomb penalty) and to three lower en- 
ergy tig levels (spanned by the three states \d xy ), \d xz ) and 
\d yz ) in which the wave-functions point away from the sur- 
rounding oxygen atoms). Following 0], the latter three t 2g 
states will henceforth be denoted by \X), \Y), and \Z). 

The super-exchange Kugel-Khomskii Hamiltonian (KK) 
l30ll depicts a coupling between orbital and spin degrees of 
freedom in transition metal 3d systems such as the t 2g titanate 
LaTi03 and the vanadate LaVC>3. As demonstrated by KK, in 
orbital systems the super-exchange Hamiltonian is anisotropic 
with a strong dependence on the orbital state. The overlap be- 
tween the transition metal orbitals amongst themselves (and 
more notably via intermediate oxygen orbitals) strongly deter- 
mines the strength (and viability if any) of the spin exchange 
interactions along the crystalline directions. In e g systems, 
the orbital component of the KK Hamiltonian exhibits, on a 
cubic lattice of size L x L x L, a discrete [Z 2 ] 3L (d = 2) 
gauge-like symmetry correspondingto reflections of spins in 
entire planes about certain axis 01, HI. As emphasized in 01, 
the three level t 2g version of the KK Hamiltonian, which is 
relevant to the titanates, prohibits direct hopping of electrons 
of a certain orbital flavor along one lattice direction. The key 
feature of the super-exchange Hamiltonian is that it prohibits 
hopping via intermediate oxygen p orbitals between any two 
electronic states of orbital flavor a (a = X, Y , or Z) along 
the a axis of the cubic lattice (see Fig[l]). 

In these materials the KK Hamiltonian l30ll can be written 
as H = H x + H Y + H z with 

H * = J Yl 4,p,**n,v4™ c M>°' (13) 

(ij)£ct /3,7/q; crrj 

In this Hamiltonian, (ij) £ a is a nearest neighbor bond along 
the cubic lattice a axis, c\^ a creates an electron at the cubic 
lattice site i in a (3 orbital state of spin a. 

As is evident from Ea.(IT3b. a uniform rotation of all spins, 
whose electronic orbital state is \a), in any given plane (P) 
orthogonal to the a axis 

cL = E^ >d U (14) 

?7 

with a, r] the internal spin polarization directions, leaves the 
KK Hamiltonian invariant. 

This leads to a conservation of the net spin of the elec- 
trons of orbital flavor \a) in any plane orthogonal to the cubic 
a axis. This manifests itself via a continuous SU(2) pla- 
nar gauge group. The continuous planar rotation Op ;a = 
[exp(z5g -6%) /h] is generated by Sf> = £ ieP 5f , the sum of 

all the spins S t,a in the orbital state a in any plane P orthog- 
onal to the direction a, is a symmetry operation (see Fig[U. 
Our theorem immediately prohibits a finite magnetization in 
these systems {S)= (Corollary III)- coinciding with the 
central result of 0]. We note that this symmetry does not 



prohibit the much more robust nematic spin ordering. For in- 
stance, in the presence of orbital ordering in the \a) state, the 
order parameter P = (SV ■ ^V+e^ ) along the 77 = x,y, or z axes 
with r] 7^ a need not vanish by our theorem. In real systems, 
the ideal mathematical model of KK with its perfect gauge- 
like symmetry and Elitzur theorem forbidden magnetization 
is likely corrupted by small effects (e.g. the hopping ampli- 
tude t± along the a axis does not identically vanish). Indeed 
finite temperature magnetization is seen in many of these com- 
pounds. Generically, in systems where the gauge-like symme- 
tries inhibiting order are lightly lifted, we may naturally antic- 
ipate the now allowed ordering to be fragile and to exhibit a 
low transition temperature. 





x( y^^J 



FIG. 1 : (color online) The anisotropic hopping amplitudes leading to 
the Kugel-Khomskii (KK) Hamiltonian fEa. J13t ). Similar to U, the 
four lobed states denote the 3d orbitals of a transition metal while 
the intermediate small p orbitals are oxygen orbital through which 
the super-exchange process occurs. The dark and bright shades de- 
note positive and negative regions of the orbital wave-function. Due 
to orthogonality with the intermediate oxygen p states, in any orbital 
state I a) (e.g. \Z) = \d xy ) above), hopping is forbidden between 
two states separated along the cubic a (Z above) axis. The ensu- 
ing super-exchange (KK) Hamiltonian exhibits a two-dimensional 
577(2) gauge symmetry corresponding to a uniform rotation of all 
spins whose orbital state is \a) in any plane orthogonal to the cubic 
direction a. In the above all spins whose orbital states are \Z) have 
been rotated in an arbitrary plane orthogonal to the cubic Z axis. 

b) Orbitals- Orbitals in 3d systems interact amongst them- 
selves via Jahn-Teller distortions which for much the same 
reasons as in the KK Hamiltonian are anisotropic. The canon- 
ical prototype of all orbital interactions is the orbital compass 
model 13111 . To streamline the essential physics in its sim- 
plest form, we examine the two-dimensional version of the 
orbital compass model. A complete mathematical analysis is 
detailed elsewhere [10]. In the two-dimensional orbital com- 
pass model at each site r of an L x L square lattice there is a 
S=l/2 operator denoted, throughout for XY spins, by a lower 
case Sff= !<Tf (the upper case S will be reserved for three 
component spins). The planar orbital compass model Hamil- 
tonian is: 
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In the three-dimensional model, additional nearest neigh- 
bor o z o z bonds appear parallel to the z-axis. We will come to 
this in Eq.i2Tt later on. To tabulate the symmetries in this sys- 
tem, we define an operator on an arbitrary line P (of intercept 
y) Op ;y = rirGP a f? W ^ m a similar definition for Op ;x . Up 
to global multiplicative phase factors, Op- y is a rotation by 7r 
about the x axis of all sites r in the line P whose x component 
is r x = x. Similarly, the operator Op- x = exp[i(7r/2)crp] 
with op = XlreP a f' This is very much similar to the op- 
erators which appeared in our discussion above of the planar 
gauge-like symmetries of the KK model, yet for the orbitals 
here, these symmetries are discrete (the angle of rotation Op 
is quantized). From the products Op^ z a^ ,x,y Op- z = —crp y , 

it is clear that Op ;z is a symmetry of H . Similarly, Op- x 
inverts the y-component of all spins on the line of intercept x 
while leaving a x untouched. In arbitrary spin size (S) versions 
of this model, the string operators above containing spin- 1/2 
Pauli matrices may be replaced by the corresponding spin-S 
7r-rotations of all spins in the line P . These "string" oper- 
ators spanning spin operators on entire lines commute with 
the Hamiltonian, [H, Op- a ] = 0. As in the case of spins in 
orbital ti g systems, rotations of individual lower-dimensional 
"planes" (lines) about an axis orthogonal to them leave the 
system invariant. As a consequence of these symmetries, each 
state is, at least, 0(2 2L ) degenerate. 

As there exist one-dimensional gauge-like symmetries 
Op-,x,y in this system, by corollary (ii) the magnetization van- 
ishes identically and only quantities invariant under Op ;XyV 
may attain a finite expectation value. Indeed, it has been es- 
tablished that a nematic-type order respecting the gauge-like 
symmetries exists in related classical extensions of this sys- 
tem OOlEi. More sophisticated order parameters (T) (in- 
variant under {O^}) and their composites may be triv- 
ially constructed 11011 . In these, (T) = ([{ the num- 
ber of {cnS} on any given horizontal row and the number of 
{crH,} on any given vertical column are both even. Due to the 
anti-commuting nature of and a^, any specific operator Ti 
commutes with all 2L inversion operators {Op- x , y }, with the 
Hamiltonian H and with any other operator iflOll. Ob- 
viously, the largest and most physically relevant order param- 
eter is the minimal gauge invariant dimer like nematic order 
parameters {SfSf + ^ a ) (with a = x,y), and their linear com- 
binations. 

For the three-dimensional "120° model" encapsulating the 
Jahn-Teller distortions and the orbital only component of the 
orbital-spin KK Hamiltonian in "e^" orbital systems (such as 
e.g. the colossal magneto-resistive manganate LaMn03, the 
cuprate KCUF3, the nickelate LiNi02, Rl^CrCU, and many 
other systems) novel discrete planar gauge-like symmetries 
and a finite on-site magnetization were found The 
three-dimensional "120° model" model possesses an exact 
discrete (d = 2) \Z^ L gauge-like symmetry (by sequential 
planar Rubick's cube like reflections about internal spin di- 
rections; See Fig.l in Refs. for an illustration). Unlike 
the case for the spins in the KK Hamiltonian for orbital t^g 
systems, there are no exact nor emergent continuous gauge- 
like symmetries. In accord with our postulates, the three- 



dimensional 120° model can (and indeed does) have a finite 
magnetization |0], Ol- 

C. Exact Gauge-Like Symmetries: Superconducting Arrays 

A related system having only four-spin interactions (i.e. 
with absolutely no two spin interactions), which emulates (p+ 
ip) superconducting grains (such as those of e.g. S^RuC^) on 
a square grid, was recently investigated by Xu and Moore @], 
|5|] . The Hamiltonian of this system is 

H = -K^a z a z a z a z -h^a*. (16) 

□ r 

Here, the four spin product is the product of all spins common 
to a given plaquette □. Note that the spins reside on the ver- 
tices on the plaquette (not on its bonds as gauge fields). It has 
been established via duality mappings iflOll. that this model 
of superconducting arrays is identical to the two-dimensional 
orbital compass model discussed in IIIBb). Not too surpris- 
ingly, the systems investigated by Xu and Moore have clear 
one-dimensional gauge-like symmetries similar to those of the 
two-dimensional orbital compass discussed in IIIBb). Rather 
explicitly, 

Op = n cj% (17) 

feP 

with P any horizontal row or vertical column, is a symme- 
try of Ea.dT^l. These symmetries are trivially related to the 
gauge-like symmetries of the planar orbital compass model 
upon invoking the duality mappings of lHoll . 

D. Other systems with Gauge-Like Symmetries 

Many other pioneering works unveiled plenty of new 
gauge-like symmetries in other systems. In B13I1 . it was 
shown that ring exchange Bose metals, in the absence of 
nearest neighbor Boson hopping, exhibit much the same one- 
dimensional gauge-like symmetries discussed above. Similar 
sliding symmetries of Hamiltonians (actions) invariant under 
arbitrary continuous deformations of a certain field along a 
transverse direction, 

<l>(x,y)^<l>(x,y) + f(y), (18) 

appear in many systems. As much throughout this work, in 
these systems there exists no rigidity to shear deformations 
along one (or more) direction(s) (y in Ea.(IT8l). Many such 
systems were discovered and discussed at length in rather 
novel works on Quantum Hall liquid crystalline phases Q, 
fl5ll . a number of models of lipid bilayers with intercalated 
DNA strands lH^I . and "Luttinger Liquid Arrays" or "sliding 
Luttinger liquids" iTTvIl which, as suggested in various works, 
may be intimately related to "stripes" 1321 . 133B . The conse- 
quences of our theorem (and Corollaries II, III in particular) 
are immediate for these systems, (<j)) = 0. 
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E. Exact Gauge-Like Symmetries by Inverted Dimensional 
Reductions 

Many spin systems can be systematically generated to have 
the gauge-like symmetries of orbital systems (as well as many 
others) |6]. We found a very simple algorithm doing so which 
very lucidly invokes dimensional reduction in the reverse or- 
der. The basic idea underlying our scheme is that a high- 
dimensional system containing decoupled chains or planes 
is, of course, one-dimensional and possesses trivial one- or 
two-dimensional gauge-like symmetries. To illustrate, we 
first consider a one-dimensional XY spin chain having nearest 
neighbor interactions, 

H = jY j ^-s 3 , (19) 

(ij) 

with i and j nearest neighbor lattice sites and s= (s x ,s y ) an 
XY spin. This translationally invariant Hamiltonian may, of 
course, be expressed in Fourier space, 

H =^^s a (k)v a/ 3(k)sp(-k), (20) 

k 

where a and f3 denote the two internal spin polarization di- 
rections, and the matrix v a p = (2 J5 a p cos k). Let us now 
make this one-dimensional system higher-dimensional by re- 
placing the scalar quantity k by more complicated quantities. 
If we replace the scalar momentum k of the one-dimensional 
spin chain by a two-dimensional k which appears in v(k) 
only via [cos(fc • a)] with a a real vector, then the two- 
dimensional system will be a trivial realization of decoupled 
chains (interacting along the direction a) and trivially pos- 
sessing one-dimensional gauge symmetries. We may now 
replace the c-number coefficients a$ by operators to obtain 
non-trivial systems. Some of the simplest examples are af- 
forded when {cii} are chosen to be projection operators. For 
instance, by replacing the scalar k of the one-dimensional 
spin chain by a simple two-dimensional diagonal matrix, we 
obtain v(k) = 2Jdg(cos k x , cos k y ). (Here and throughout, 
dg(ai, a2, ad) refers to a d-dimensional diagonal whose 
diagonal elements are {ai, a2, ad})- In the v(k) above, 
we replicate the single momentum coordinate to generate a 
higher-dimensional interaction. Insertion of this matrix into 
Ea.(l20l. leads to the two-dimensional orbital compass model 
of Bb) having [Z 2 ] 2L gauge-like symmetries of "en-block" ro- 
tations by 7r of all spins belonging to a certain row/column. 

Similarly, by setting v(k) = 2Jdg(cosk x ,cosk y ,cosk z ) 
in Eq.(l20l» for a three component spin system S = 
1 Sy-) S xz ), we obtain the three-dimensional orbital com- 
pass model (on a cubic lattice of size L x L x L) 

H = J^2(SfSf +ix + SfS^ + S*S* +i J. (21) 

This system has an exact [Z2] 3Ij2 symmetry for classical spins 
(along each chain parallel to the cubic a axis, we may reflect 
S a — > — S a yet keep all other spin polarizations unchanged, 



SpjL a — > Sp-L a ) as well as a lower exact [Z2] 3L gauge-like 
symmetry (forming a subset of the larger [Z2] 3L symmetry 
present for classical spins) when the spins are quantum: in 
this case we may rotate all spins in a plane orthogonal to the 
cubic lattice direction a by tt about the internal S a quantiza- 
tion axis. 

In Ea.(l20l» for an XY spin system, we may further gen- 
erate interactions along arbitrary tilted rays by considering 
v = 2Jcos(k x + Xk y as) with as the Pauli matrix. An in- 
sertion this kernel into Ea.(l20l> leads to real space interactions 
of the x components of the spins along chains parallel to the 
(1, A) direction and of an exchange interaction amongst the y 
components of the spins along the (1, — A) direction. By tun- 
ing the A — > 0, the two "clapping" (1,±A) directions fold 
back and degenerate into the original one-dimensional XY 
spin chain. We find that for all A / 0, the classical sym- 
metries are much the same as for the two-dimensional orbital 
compass model: reflecting all s x — > — s x while leaving s y un- 
touched along any chain parallel to (1, A). Thus, according to 
our theorem and corollary (II), no finite on-site magnetization 
is possible in classical XY systems with the above Hamilto- 
nian. Nematic like order is possible and was indeed found 
0] for A = d=l. For A ^ ±1, there is, however, no canonical 
transformation that may be applied to the quantum spin model 
to encapsulate these symmetries. 

If we start with a nearest neighbor three component 
(Heisenberg) spin model on the square lattice and re- 
place the usual isotropic interaction kernel 2Jdg(cosk x + 
cos ky, cos k x + cos k y ) by v = dg(cos k x + cos k y , cos k x + 
cos ky , cos k z ) to generate a spin system residing on a cubic 
lattice, we obtain a new system with many exact gauge-like 
symmetries. In classical systems, for all chains parallel to 
the z-axis we may take S z — > —S z and for all planes par- 
allel to the XY plane we may perform a uniform continuous 
rotation about the z axis- a [U(1)] L x symmetry. If the 
spins are quantum, only the canonical [£7(1)] L symmetries are 
present. Our theorem precludes finite on-site order for both 
classical and quantum spins on a cubic lattice with the above 
spin Hamiltonian. 

We may similarly start with many other lower-dimensional 
systems to generate higher-dimensional systems by these op- 
erations. In some instances, these transformations lead to sys- 
tems of intense physical interest (e.g. orbital compass mod- 
els employed for elucidating the properties 3d transition metal 
compounds at large |31], and crystal field split t^g systems in 
particular). For each of these systems, our theorem and its 
corollaries can be applied to often inhibit the possibility of a 
finite on-site magnetization. It should be emphasized that al- 
though these models are generated by artificially "lifting" a 
low-dimensional momentum, the resulting Hamiltonians are, 
obviously, not lower dimensional. We are only considering 
symmetry properties here and not the thermodynamics. In the 
current context, "dimensional reduction" may be understood 
from a technical point of view as an algorithm for yielding 
gauge-like symmetries in a natural fashion. Thermodynamic 
dimensional reduction of these models occurs only in the large 
n (spherical) limit of classical renditions of all of the above 
mentioned spin systems (e.g. the two-dimensional compass 
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model is identical to an XY chain). We will briefly com- 
ment to other special properties of large n systems when dis- 
cussing new emergent gauge-like symmetries that we found 
in Brazovskii-type systems relevant, amongst other things, to 
liquid crystal systems, theories of structural glasses, electronic 
systems favoring phase separation that is inhibited by long 
range Coulomb interactions, and to many other arenas. 

By reverse engineering various given gauge-like lattice 
symmetries to find corresponding Hamiltonians, we found ad- 
ditional models on different lattices having gauge-like sym- 
metries when polarization (spin space) dependent spin-spin 
interactions were slaved to the bond directional on the lattice 
(external space) to produce many different real space inter- 
action kernels v a p(k) (or their direct real space transforms 



F. Emergent gauge-like symmetries in spin, orbital, and other 

systems 

We now examine situations wherein a gauge-like symmetry 
emerges within the low energy sector. In some instances, it is 
present in leading order fluctuations about the classical ground 
state (e.g. 1/5 spin wave corrections). From corollary (IV), 
we know that in the presence of a gap, expectation values of 
any symmetry non-invariant quantities must vanish. In gen- 
eral systems having various emergent gauge-like symmetries, 
non-gauge invariant quantities may often vanish. I. Quantum 
systems: a) Geometrically Frustrated Spin systems: In many 
spin systems (e.g. the kagome, pyrochlore, and checkerboard 
magnets), the availability of many spins vis a vis the number 
of building blocks on which they reside allows the prolifer- 
ation of many low energy states (3. Pyrochlore magnets 
such as the S= 3/2 spinel ZnG^C^ have remarkable magnetic 
behavior (such as e.g. very clearly detectable spin waves link- 
ing degenerate states). The checkerboard lattice ll35ll may be 
viewed as a two-dimensional projection of the three dimen- 
sional pyrochlore lattice composed of corner sharing tetrahe- 
dra on the plane. The system is composed of nearest neighbor 
bonds on the square lattice augmented by additional diagonal, 
next nearest neighbor, interactions along all diagonals of one 
plaquette sub-lattice (hence the name "checkerboard"). Both 
the pyrochlore and checkerboard magnets has an exponential 
in volume classical ground state degeneracy as in other sys- 
tems of emergent local (or d = 0) gauge symmetry. Lead- 
ing order (1/5) corrections were computed finding the effec- 
tive spin- wave Hamiltonian for fluctuations about the ground 
states (35l l^]. Much like the orbital systems discussed in 
Section (IIIBb), the resulting Hamiltonian possesses one di- 
mensional gauge like symmetries. In l37ll . an effective low 
energy Hamiltonian was derived for these quantum systems 
without the need for spin- wave theory. The resulting effective 
Hamiltonian looks much the same as the KK Hamiltonian for 
orbital systems discussed earlier and displays much the same 
discrete gauge-like symmetries seen in the KK relevant to or- 
dering in e g materials 101,101 • Gauge-like theories for geo- 
metrically frustrated magnets with string operators generating 
freely propagating "spinons" are also found elsewhere, e.g. 



[38] for the pyrochlore magnet. The large 5 ground states of 
kagome antiferromagnets may be linked to each other via all 
possible contortions of a two-dimensional membrane. The ef- 
fective low energy Hamiltonian of the quantum Kagome mag- 
net is, once again, similar to the KK Hamiltonian. 11391 

b) Spin systems with ring exchange: In a recent work 
emergent one dimensional gauge like symmetries were found 
in a quasi-exactly solvable square lattice spin system of a 
Heisenberg antiferromagnet augmented by ring exchange in- 
teractions of comparable strength. Here, the presence of one- 
dimensional gauge like symmetries went hand in hand with 
fractionalization: spinons propagating along the diagonal di- 
rections incurred no energy penalty while spin motion in other 
directions led to a confining potential. In an upcoming work 
lEoll . we identify gauge-invariant order parameters in this sys- 
tem and further suggest the existence of a finite temperature 
"order out disorder" transition suggesting natural finite tem- 
perature fractionalization in these two-dimensional systems. 

II. Classical spin systems: a) Within the classical (large 
S) limit of d=3 orbital compass of Eq.(l2Tl (relevant to such 
compounds as YVO3 LaVC>3, and LaTiOs), we find an emer- 
gent continuous rotational symmetry of all the spins in an en- 
tire plane about a spin axis orthogonal to the plane direction. 
For instance, taking any uniform spin state (automatically a 
ground state of Eq.(l2Tl) and rotating all of the spins in a given 
plane orthogonal to the z-axis (i.e. all spins on sites r satisfy- 
ing r z = z) by an arbitrary angle about the S z spin polar- 
ization axis leads to another ground state. This can be further 
mutated by similar operations in other planes (parallel and or- 
thogonal). Compounding all of the symmetries together, this 
corresponds to an emergent [U(1)] 3L degeneracy. Indeed, it 
was found (|9]) that the system has nematic like order invari- 
ant under this continuous planar gauge-like symmetry, b) The 
Fully Frustrated XY model on a honeycomb lattice: within 
the ground state of such an XY system having half a fluxon 
threading each hexagonal "plaquette" EUl . each "plaquette" 
has exactly two bonds having zero phase gradient (similar to 
the requirement of one singlet per square plaquette in H3o 
with all other bonds having a phase gradients of 7r/4. Con- 
necting the zero gradient bonds to each other leads to non- 
intersecting stripes threading the system. There are several 
ways of generating such stripe configurations, all of which are 
inter-related by an emergent one dimensional like gauge sym- 
metry, c) An Ising spin systems on a s quar e lattice with com- 
peting next and next nearest neighbor 14211 interactions, hav- 
ing a shell of minimizing modes in Fourier space, was found 
to possess an one-dimensional emergent gauge-like symmetry 
with a diagonal stripe like structure within its ground state, d) 
General two spin Fourier kernels, such as v(k) = (k 2 — q 2 ) 2 
in Eq. (l20l) . with a (d — 1) dimensional shell of minimizing 
modes in Fourier space are relevant amongst other things to 
the Brazovskii models l43ll for liquid crystals, block copoly- 
mers, frustrated electronic systems, and structural glass tran- 
sitions |3 We may rigorously establish, in the large 
n (spherical) limit, that the ground state entropy scales as 
the surface area of the system, 5 oc (qN)^-^^ |46]. The 

ground state possess an SO(N < ^ d ~ 1 ^ d ) symmetry of which a 

N (d-i)/d 

gauge-like [U(l)] 2 is a subgroup. The Brazovskii sys- 
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tem is endowed with a continuous symmetry harboring a one 
dimensional gauge-like character. Similar to our corollary for 
gapped quantum systems, in the large n Brazovskii system 
the on-site magnetization is zero and, in fact, the system pos- 
sesses a zero temperature critical point. Indeed in all trans- 
lationally invariant systems with a kernel v(k) analytic about 
its minima and in which the minimizing Fourier modes form 
a (d — 1) dimensional manifold (as in rotationally invariant 
incommensurate systems) or a (d — 2) dimensional manifold 
corresponding to continuous one or two-dimensional gauge 
symmetries. In large n spin systems, any permutation of the 
values of the kernel v(k) — > v(Pk) (with P any permutation 
of the TV momenta) leaves the partition function identically 
invariant B47I1. This allows an exact reduction to decoupled 
chains at low temperatures where gauge-like symmetries triv- 
ially exist. Such symmetries were also shown to rigorously 
exist in large n renditions of vectorial spins when subjected 
to uniform non-Abelian gauge backgrounds. EUl [Within the 
spherical model, the inverse critical temperature is given by 



k B T c J (2n) d v (k)- v (g) 

which diverges for (d— 1) or (d — 2) dimensional manifolds of 
modes q which minimize an analytic kernel v(k).] These sym- 
metries lead to "zero-energy domain walls" which go hand in 
hand with fractionalization in the quantum models of fl2ll . 
These high degeneracies do not generically survive for finite 
n system (where, for example, XY spin systems with the 
same kernel have an emergent discrete Zm symmetry (with 
M ~ (gL)f _1 )). Nevertheless, in some special instances they 
do, e.g. O. 

IV. DISCUSSION 

The degeneracy found in all (or only ground) states in sys- 
tems having exact (or, respectively, only emergent) gauge-like 
symmetries suggests deep non-trivialities in the ordering tran- 
sition that occurs, if any, at low temperatures. 

As well known, in the presence of ground state degenera- 
cies, entropic fluctuations often drive the system to order in 
the vicinity of "softer" states having more phase space for 
low energy fluctuations about them |3- This mechanism 
is often termed "Order out of Disorder". The origin of this 
name are the entropic fluctuations ("disorder") which drive the 



low temperature ordering transition. The presence of gauge- 
like symmetries makes many standard methods of analysis 
much harder or inexact. Nevertheless, finite temperature or- 
bital ordering (via an "order out of disorder" mechanism) was 
rigorously established in @], |0] for two classical systems 
having gauge-like symmetries. Finite temperature magneti- 
zation was proved in the 120 degree model which harbors 
an exact [Z2] 3L symmetry. Nematic order was also shown 
to rigorously exist in the classical three dimensional orbital 
compass model which has an exact \Z^ L symmetry and an 
emergent [U(1)] 3L symmetry. The techniques presented in 
|0], 131 may be straightforwardly applied to a host of other 
systems possessing gauge-like symmetries, e.g. the effec- 
tive low energy Hamiltonians of H37ll . ll39ll for the pyrochlore 
and kagome magnets. The natural candidates for order pa- 
rameters are quantities invariant under all low dimensional 
gauge-like symmetries (exact or emergent). With this guid- 
ing principle in hand, a natural order parameter is found for 
two-dimensional quantum spin systems lfl2il exhibiting finite 
temperature fractionalization. H4QB 

In summary, we showed that the generalization of Elitzur's 
theorem to d-dimensional gauge groups (intermediate be- 
tween local and global gauge symmetries) provides a formal 
meaning for the notion of dimensional reduction. Our the- 
orem proves that for quantities which are non-invariant un- 
der the d-dimensional gauge transformations, the effective di- 
mension of /^-dimensional theory under consideration is re- 
duced from D to d. Elitzur's theorem is then a consequence 
of the absence of phase transitions at any finite temperature 
in a zero-dimensional system (see our corollary I). Moreover, 
this result also implies that a one-dimensional gauge symme- 
try cannot be broken if the interactions of the theory have a 
finite range (corollary II). The same holds true for continuous 
two-dimensional gauge symmetries (corollary III). We further 
described different theories of current interest which possess 
exact gauge-like symmetries for which corollaries II and III 
directly apply. Additional systems are known to have emer- 
gent gauge like symmetries and we illustrated (corollary IV) 
how Elitzur's theorem extends even to such systems in the 
presence of gaps. The wealth of systems of current interest 
for which our theorem and its four corollaries directly apply 
illustrates the power and generality of our result. 
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k 



with (3 the inverse temperature. As is evident from Eas. J23l24l . 
for any of the AM permutations P of the spectra {v(k)} — > 
{v(Pk)}, the partition function Z is invariant. 



with fx determined from 

p= I (24) 
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